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ABSTRACT 

Carbon-rich evolved stars from the asymptotic giant branch to the planetary nebula phase 
are characterized by a rich and complex carbon chemistry in their circumstellar envelopes. A 
peculiar object is the preplanetary nebula SMP LMC 11, whose Spitzer-IRS spectrum shows 
remarkable and diverse molecular absorption bands. To study how the molecular composition in 
this object compares to our current understanding of circumstellar carbon chemistry, we modeled 
this molecular absorption. We find high abundances for a number of molecules, perhaps most 
notably benzene. We also confirm the presence of propyne (CH3C2H) in this spectrum. Of all 
the cyanopolyynes, only HC3N is evident; we can detect at best a marginal presence of HCN. 
From comparisons to various chemical models, we can conclude that SMP LMC 11 must have an 
unusual circumstellar environment (a torus rather than an outfiow). 



Subject headings: astrochemistry- 
stars: individual (SMP LMC 11) 

Introduction 



As stars with initial masses between 0.8 and 8- 
9 Mq approach the end of their lives, they reach 
the asymptotic giant branch (AGB) stage. The 
AGB is characterized by alternate hydrogen and 
helium shell burning, dredge-up events and ex- 
treme mass-los s rates (up to 10~^ Mq year~^, 
Iben fc Ren^ [l983l. This high mass loss causes 
the star to evolve further, leaving the AGB. When 
a star leaves the AGB it may become what is 
known as a preplanetary nebula (pPN) during 
a relatively short-lived ( lasting '^lO'^-lO'* years. 



-circumstellar matter — stars: AGB and post-AGB — stars: carbon- 



Vassiliadis Sz W00J1994 transitionary period be- 



fore becoming a planetary nebula (PN). 

The material lost by a star during its time on 
the AGB goes into the circumstellar environment 
(CSE) before it is dispersed into the interstellar 
medium (ISM). The CSE is a relatively cool re- 
gion (with temperatures lower than the effective 
temperatures of AGB stars, where Toff ~ 3000 
K), which allows the formation of molecules and 
dust beginning on the AGB and continuing into 
the (p)PN stages. 

One of the first and most stable molecules to 
form in the CSE is CO. As a result, the rela- 
tive amounts of carbon and oxygen in the CSE 
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largely determines future chemistry. Stars begin 
their lives with more oxygen than carbon (C/0 < 
1), but depending on the initial stellar mass and 
metallicity, an AGB star may undergo sufficient 
dredge-up events and become carbon-rich (C/0 > 
1), resulting in what is known as a carbon star (or 
a carbon-rich star). 

The chemistry of carbon stars can result in an 
assortment of molecules due to the ability of car- 
bon to form a variety of chemical bonds. For ex- 
ample, more than 60 molecules have been detected 
in the CSE of the pr ototypical carbon-rich AGB 



star IRC-K10216 (e.g. ICernicharo et al.lll996f ) 



The carbon chemistry in the CSE of AGB stars 
and (p)PNe is also of particular interest because a 
group of molecules known as polycyclic aromatic 
hydrocarbons (PAHs) are thought to form here 
(Latter 199lh . These molecules are ubiquitous in 
the universe with u p to 10-15% of cosm ic carbon 
contained in PAHs (|Snow fc wit^ll995[ ). 

In the current models of PAH formation, ben- 
zene (CfjHfi) formation is considered a bottleneck 



step ( Frenklach fc Feigelsonll989l : ICherchneff et al. 



19921 ^. 

Not much is known about the formation of ben- 
zene in the CSEs of evolved stars, although there 
are chemical models describing its formation in 
the inner shock ed regions of the CSE surround- 



ing AGB stars (jCherchneff et al.lll992l l as well as 



throug h photochemical reactio ns in the CSE of 
pPNe (|Woods et al.ll2002ll2003l) . 

While PAHs are ubiquitous in the universe (in- 
cluding carbon-rich PNe), benzene is not often 
found in evolved stars; it has b een found in just 
two ob jects thus far: CRL 618 ( Cernicharo et al 
200 lal) and SMP LMC 11 (|Bernard-Salas et al 



2006L hereafter Paper I). This suggests that it is 
either difficult to form or that it reacts quickly 
once formed. 

In this paper, we will discuss the latter object, 
SMP LMC 11, which is a carbon-rich pPN in the 
Large Magellanic Cloud. It is d escribed as a low- 



excitation pPN (jSanduleak et a l j 119781: iMorgan 



[1984) and has a bipolar outflow (|Shaw et alll2006[ ) 



It also shows a rather high expansion velocity (122 
km/s, iDopita et al. 1988) with multiple velocity 
components. 

Here we present a detailed analysis of the 
molecular absorption bands in the mid-infrared 



spectrum of SMP LMC 11 (first presented in Pa- 
per I) and show that our results are inconsistent 
with current models for the chemistry in evolved 
carbon-rich CSEs. 

We begin this paper with a description of the 
observations and data reduction in Section [2l then 
we describe the dust continuum in Section [3l We 
follow this with an inventory of the molecular 
bands in the spectrum and the method we use to 
model these bands in Section S] We then present 
our results from our model fits in Section [5] Next, 
we discuss the implications of our results for the 
evolutionary status, chemical evolution and geom- 
etry of SMP LMC 11 in Section El Finally, we 
present our conclusions in Section [71 

2. Observations and Data Reduction 

SMP LMC 11 was observed with t he Spitzer 
Space Telescope (IWerner et al.l 2004) Infrared 
Spectrog raph (IRS. lHouck et al.ll2004l) as part of 
the GTO program on June 6, 2005 (program ID 
103, AOR key 4947712). Here we present a new 
reduction of the spectrum with the latest calibra- 
tion files (pipeline version S18.18). We obtained 
the basic calibrated data (BCD) files for SMP 
LMC 11 and processed the data for the short high 
(SH, R = 600, A = 9.9-19.6 short low (SL, 

R = 60-127, A = 5.2-14.5 /zm) and long low (LL, 
R = 57-126, A = 14.0-38.0 /im) modes. 

We cleaned the data using irsclean with the 
campaign rogu e pixel mask and extracted it in 



SMART v8.2.1 ( Higdon et al.l 120041 ): we extracted 



the SH data using full aperture extraction, and 
the SL and LL data with the manual o ptimal ex- 
traction mode ( Lebouteiller et al. 2010() . Next we 



defringed the LL mode and trimmed the edges of 
the orders for all modules to remove edge effects. 
We eliminated flux jumps between the orders by 
comparing the overlap regions and scaling the or- 
ders (for the SH, orders were scaled to match order 
20, for the low resolution data SL2 was scaled to 
SLl and both were scaled to the LL data), then 
we averaged the flux from the two nod positions. 
Since we were unsure of the reliability of the ini- 
tial uncertainty estimates, we instead estimated 
the uncertainties on the flux values by measuring 
the standard deviation in a featureless region of 
the SH spectrum between 16.53 and 17.44 ^m; we 
found a standard deviation of 0.0122 Jy, corre- 
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SMP LMC 1 1 with Dust Model 




15 20 25 

wavelength (ji/m) 



Fig. 1. — The Spitzer-IRS low- resolution spectrum 
(SL and LL combined) of SMP LMC 11 (black) 
and a dust model (dashed grey or dashed blue 
line in the electronic edition). The inset shows 
the 2MASS fluxes with the dust model. See the 
electronic version of the Journal for a color figure. 

sponding to a S/N ratio of 47 in this range. 

Finally, to facilitate comparison to molecu- 
lar models, we shifted the spectrum to the rest 
frame using a radial v elocity of 263.5 km s~^ 
(|Morgan fc ParkeilllQQSh and the relative motion 
of Spitzer at the time of the observations (Vlsr 
= 12.5 km/s). The full low resolution spectrum is 
shown in Fig. [TJ 



3. The dust continuum 

Fig. [T] shows the Spitzer-IRS continuum as well 
as the 2MA SS (ISkrutskie et al.. 2006) H, J and K 
band fluxes ( Cutri et al.l 20031) From those data, 
it seems clear that there is no appreciable stel- 
lar continuum that could contribute to the mid-IR 
emission. We thus conclude that the fiux in the 
entire IRS spectrum is dominated by dust emis- 
sion. 

The region beyond 17 /im has a few apparent 
features near 19 and 21 fim in the LL data. How- 
ever, as these features are absent from a quick re- 
duction of the long high (LH, R = 600, A = 18.7- 
37.2 fim) data, we consider these to be artifacts. 



1.2 1.4 1.6 1.8 2.0 2.2 
wavelength (firr)) 



and an ar- 
""''^ of 25 



Thus, the dust emission does not show clear spec- 
tral features, and thus is presumably caused by 
amorphous carbon dust. 

Since it is precisely this dust emission that is 
then absorbed by molecular bands (especially in 
the 10-17 /im region, see Fig. [ij, we require some 
idea about the properties of this dust continuum 
before we can model the molecular bands. A sin- 
gle blackbody curve cannot properly reproduce the 
overall shape of the dust continuum, suggesting 
that there is some stratification in the dust layers. 
From our point of view, we can see through the 
outer (and colder) dust layers up the point where 
the dust becomes optically thick; at that point, the 
dust temperature is T^^^. We thus approximated 
the dust emission by a weighted sum of blackbody 
spectra at temperatures between 
bitrarily chosen minimum temperature 
K in steps of 25 K. Assuming optically thin dust 
in radiative equilibrium with the hottest (optically 
thick) dust, we then determined the appropriate 
weights for each temperature bin consistent with 
a constant- velocity outfiow at a constant mass-loss 
rate (i.e. we determined the mass in each layer). 
We find that the best model is one with a maxi- 
mum dust temperature of about 425 K (see Fig.[T]). 
However, even in such a case, we overestimate the 
dust emission at the longer wavelengths; this sug- 
gests that there is less cold dust than expected for 
a constant- velocity outflow at a constant mass-loss 
rate. Nonetheless, the cold dust does not matter 
much for the analysis of the molecular bands; for 
the modeling presented in this paper, we assume 
that the molecules are absorbing 425 K blackbody 
radiation. 

We also note that our dust temperature esti- 
mate is higher than that in Paper I, but our esti- 
mate should be more reliable than their greybody 
flt. 

4. Molecular bands 

4.1. Molecular Inventory 

Superposed on the dust emission are many ab- 
sorption features due to various molecular species 
(see Fig. (2). We will restrict our discussion to 
wavelengths shorter than 17 fim, where most of 
the absorption occurs. 

Here, we describe the molecular bands in great 
detail and determine the physical conditions by 
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modeling these bands. 

Acetylene (C2H2) is often the dominant molec- 
ular absorber in the infrared spectra of carbon-rich 
pPNe. In the spectrum of SMP LMC 11, C2H2 
provides by far the strongest absorption. The Q- 
branch for the 1^5 bending mode is obvious as a 
very deep absorption feature at 13.7 ^m; addition- 
ally, the P- and R-branches of C2H2 cause much of 
the broad and deep absorption that is obvious be- 
tween 12 and 16 /im (see Fig. [T]). Such broad and 
deep C2H2 absorption is also seen in some other 
carbon-rich objects, such as the carbon star IRAS 
04496-6958 (|Speck et al.ll2006[ ). for example. In 
addition, acetylene exhibits a much weaker com- 
bination band (vi+i's) at 7.5 fim which is blended 
with other features in the spectrum of SMP LMC 
11. 

Larger acetylene chains ar e less commonly ob - 
served in the CSE of pPNe (jPonfria et al.ll201ll) . 
In the spectrum of SMP LMC 11 though, diacety- 
lene (C4H2) is another major contributer to the 
molecular absorption in the spectrum. The i^s 
bending mode of this species appears at 15.9 /im 
and the z/g + i^s combination band is clearly visi- 
ble at ~8 /im. Small amounts of the corresponding 
Vii bending mode of triacetylene (C6II2) might be 
present in the red wing of the 15.9 /im feature. 

In contrast, HCN — which often shows up in car- 
bon star spectra alongside the C2II2 band at 14 
/im — at best barely contributes to the absorption 
in this spectrum. It is thus surprising that HC3N 
is present and shows a clear and deep absorption at 
15.03 /xm from its bending mode. There is no ev- 
idence for any longer cyanopolyynes; for instance, 
there is no absorption at 15.57 /im from the vj 
band of HC5N. 

Nearly unique for pPNe (wit h one other detec- 
tion thus far in CRL 618, see ICernicharo et al 



2001al) . the spectrum also shows significant ben- 
zene absorption (as noted in Paper I). It has its 
deepest absorption at 14.85 /im, where the 1/4 
bending mode absorbs more than 15% of the total 
continuum flux! Additionally, strong absorption 
is clearly present at the wavelengths where other 
strong benzene bands are expected: another bend- 
ing mode (h'li) at 9.6 /im and a ring stretching and 
deforming mode (1^13) at 6.72 /im. 

In the blue wing of the C4II2 band, near 15.78 
/xm, some additional absorption could be due to 



the 1^8 bending mode of propyne (CH3C2H, also 
sometimes called methylacetylene) at 15.78 /im, 
as suggested in Paper I. Furthermore, there is also 
some absorption visible in the spectrum from the 
1/4 band of CH4 at 7.7 /im as well as the 1^7 band 
of C2II4 at 10.53 /im. Finally, we could not de- 
termine the origin of the absorption feature near 
10.38 /xm, although we suspect a molecular origin 
for this feature. 

4.2. Modeling the molecular absorption 

We modeled the molecular absorption using the 
same methods that a re used to build the Spec- 
traFactory database ( Cami et al. 2010t ). These 
model calculations start from line lists detailing 
the frequencies and intensities of the individual 
molecular transitions. 

Line hsts for C2H2 (including the H^^CCH 
isotopologue), HCN, CH4 and C2H4 are taken 



from the HITRAN 2008 database (jRothman et al 
2009f ): the line lists for C4H2, CgHe, HC3N 
and CH3C2II are from the GEISA database 
( Jacquinet-Husson et al. 20081 ). 

As we could not find reliable line lists for all 
species or bands, we calculated some line lists 
from molecular constants and pgopher v 7.1.108 
(IWesternllioiol ). The GEISA line hst for ben- 
zene contains only data for the fundamental 1^4 
band. Thus, in order to model benzene absorp- 
tion at the shorter wavelengths in the SL data, we 
calculated line lists for the t^i3 and v^4 bands us- 
ing molecular constants from iDang-Nhu &: PHva 
(l989t ). Similarly, the GEISA line list for C4H2 
does not contain data for the transitions of the 
j/g + i^g combination band at ^8 /im; we thus cal- 
cu lated a line list using molecular constant s foun d 



Arie fc John'^ (Il992l) : ICuelachvih et all (|l984l ) 
Khlifi et al.l (Il995i r~ Fil ially, we calculated a 



C6H2 list using data from iMcNaughton fc Bruget 
(|l99ll) . 

From these line lists, we calculated optical 
depths assuming a population in local thermo- 
dynamic equilibrium and a Gaussian intrinsic line 
profile with a width of 10 km s~^, which is typical 
of outflows of evolved stars. We carried out ra- 
diative transfer through isothermal, plane-parallel 
molecular slabs in front of a 425 K blackbody 
background (see Section ^ and smoothed and re- 
binned the resulting models to match the resolving 
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power of the observations — 600 for the SH module 
and 90 for the SL data. 

We compared the resulting models to the obser- 
vations (normalized by a cubic spline continuum) 
and calculated xt^ the reduced statistic. How- 
ever, we note that there may be some systematic 
errors as well. For instance, it is important to re- 
alize that the current line lists for C2H2 do not 
allow us to reproduce the broad and deep ab sorp- 
tion in the 12-16 fim (see Speck et al. 20061 for a 
discussion). 

To find the best model, we calculated models at 
different temperatures ranging from 200 to 400 K 
with a step size of 25 K, and similarly at column 
densities between N = 10^^ and N — 10^^ cm^^ 
in steps of log N = 0.1 and then calculated the 
xt value for each model. For wavelength ranges 
containing absorption due to several species, we 
properly treated line overlap by summing the op- 
tical depth profiles for each contributing molecule 
prior to performing the radiative transfer calcu- 
lations when we fit several species simultaneously. 
We thus simultaneously modeled the absorption of 
the C2H2 isotopologues and HCN between 12 and 
14 /im and similarly also combined C4H2, C6H2, 
CH3C2H and HC3N in the 15-16.5 /im range. We 
fit benzene to the spectrum between 14.2 and 15 
fxm and C2H4 from 10.5 to 11 /im. Where overlap 
appeared across the modelled regions, we added 
the optical depth values determined from earlier 
fits to the new regions (e.g. the best fit opti- 
cal depths for C2H2 and HCN were added to the 
optical depth for the C4H2, C6H2, HC3N, and 
CH3C2H fits) and then performed the radiative 
transfer calculations. We expect that small errors 
may be introduced using this method, but due to 
the small difference in temperature between the 
layers, this effect should not be large. 

Using the results from fitting the SH observa- 
tions, we predicted the absorption in the SL data 
for benzene, C2H2 and C4H2. The CH4 absorp- 
tion was fit by itself between 7.5 and 7.9 /xm, but 
the optical depth profiles for the + vs band of 
C4H2 as well as the C2H2 band at 7.5 /xm as de- 
termined in the SH data were added to the total 
optical depth profile prior to calculating the radia- 
tive transfer calculations for the CH4 fits. 



5. Results 

The best fits for the SH and SL data are shown 
in Fig. [21 the corresponding parameters are listed 
in Table [1] As can also be seen from the residuals 
in Fig. [21 our models reproduce the observations 
quite well (xj w 0.84 over the entire fitting region). 
Our best fits indicate that the molecules are found 
in a range of temperatures from 250-375 K; how- 
ever, some of the temperature stratification may 
be artificial due to the large uncertainties on our 
temperature determinations. 

5.1. Benzene 

To fit the 1^4 bending mode of benzene at 14.85 
/im, we require a high column density (log N — 
17.80) consistent with the deep absorption in the 
spectrum. Our model reproduces the depth of the 
absorption feature well, but does not match the 
width and profile shape. This discrepancy is prob- 
ably due to the absence of hot bands in our molec- 
ular model. Hot bands are generally slightly off- 
set compared to the fundamental mode, and thus 
tend to broaden the absorption band. At the rel- 
atively high temperature of benzene found here, 
we certainly expect some contribution from the 
hot bands; at 300 K, only 54% of the molecules 



should be in the ground state (jKauppinen et al 
1980 h. However, hot bands are not included in 



the GEISA line list for benzene. As a result, our 
model fit is much narrower than it should be. 

Using the benzene parameters found from fit- 
ting the i>4 band, we predicted the appearance of 
the bands at 9.6 and 6.7 /im. As seen in Fig. [21 
while there are absorption features at these wave- 
lengths consistent with benzene, our predictions 
do not fit these bands particularly well. At 9.6 
/im, we can see that the predicted absorption is 
much deeper than the observed band, although 
the band shape is similar. 

The feature at 6.7 /im is more unusual: While 
there is an overall absorption feature in the spec- 
trum of SMP LMC 11 at this wavelength, the fea- 
ture is much broader than our model and shows 
what could be a small emission bump right at 
the central wavelength of the benzene absorption. 
Just as for the 14.85 /tm band, the broader ob- 
served feature could be a consequence of the ab- 
sence of hot bands in our model. Indeed, we note 
that the absorption in the spectrum appears to be 
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Table 1 

Temperatures and column densities for the best model fits to our data as well as 

NOMINAL la uncertainties. 



Molecule 


log N 


T(K) 


C2H4 




350 tTo 


C2H2 

H13CCH 

HCN 


18.10™ 
16 90+°-^° 
16.50_i.5o 


375 til 


CgHg 




350 


C4H2 

CH3C2H 

HC3N 


i7.ool°:J° 
i6.4oi;j:J° 


325 tlf' 


CH4 




250 



SMP LMC 1 1 Spectrum with Best Fit Models 




6 8 10 12 14 16 

wavelength (/.tm) 



Fig. 2. — The normalized SL and SH spectra for SMP LMC 11 from 6-16.6 //m with our best fit models 
and residuals. The combined best fit is shown offset above the spectrum in dark grey (red in the electronic 
edition). The individual best fits for each molecule are offset above this in light grey (blue). The residuals 
are shown offset below the spectrum and the dotted dashed lines indicate the error ranges. See the electronic 
version of the Journal for a color figure. 
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approximately twice as broad as in the model, with 
a similar red degraded wing (see Fig.[3|). Thus, we 
consider benzene the carrier for the overall absorp- 
tion feature at 6.7 /im. The nature of the small 
emission bump is not clear. If this is a real emis- 
sion feature, it seems unlikely that it would be 
due to benzene, since we do not see any similar 
emission at the other benzene absorption bands. 

5.2. Acetylene Chains 

The C2H2 absorption at 13.7 /im is the deepest 
and broadest absorption feature in the spectrum. 
Accordingly, we find the highest column density 
for this feature (log N = 18.10) of the molecu- 
lar bands we fit. Since we ignored the deep and 
broad absorption between 12 and 16 /xm which is 
at least partly due to C2H2, this is clearly a lower 
limit of the true column density. In addition to 
the main isotopologue, we detect a fairly high col- 
umn density of Hi^CCH (log N = 16.90) yielding 
a N12C2H2/NH13CCH ratio of 16 and a ^^C/^^C ra- 
tio of 33. This is in the low end of the range for 
carbon -rich objects as determined bv I Milam et al 



(|2nn9[ ). who found i^C/i^C ratios between 25-90. 
However, since we underestimate the column den- 
sity of the main isotopologue, this is most cer- 
tainly an absolute lower limit to the ^^C/^'^C ra- 
tio. Similarly, from the column densities of C2H2 
and C4II2, we find a ratio of C2H2/C4II2 of 10, 
which would again be a lower limit. 

We find that the identification of the vn band 
of C6H2 in this spectrum is uncertain. Using rea- 
sonable defringing techniques, the apparent ab- 
sorption feature near 16.1 /xm disappears (see 
Fig. HI). Further, when th e molecular constants 



from iMcNaughton fc Brugeti (|1991[ ) are used to 
calculate a model line list, we find that the cal- 
culated band center appears 0.02 fxm longer of the 
apparent absorption in the spectrum. A compar- 
ison between our calculated model and the pre- 
yiously identified C6II2 absorption in CRL 618 
(jCernicharo et al.l l2001al ) in Fig. |4] indicates that 
our calculated model is at least correct in peak po- 
sition, so we consider it possible that the apparent 
C6H2 absorption is not real. 

If we assume that the peak near 16.07 fim is due 
to C6II2 and shift our model of this band accord- 
ingly, the best fit is obtained at log N = 15.60j;Q gQ, 
which does not reproduce the shape of the band. 
In order to fit the band shape, a log N of 16.30 is 



required. However, since we have to assume that 
the position of our band as well as our continuum 
choice is incorrect to obtain this result, it seems 
unlikely that this feature is actually C6H2. 

5.3. Cyanopolyynes 

Although there is no obvious absorption from 
HCN at 14 /Ltm, our modeling attempts suggest 
that the overall fit is slightly better if HCN is 
included; in fact, we find a similar (low) column 
density for HCN as we do for the much stronger 
HC3N. To verify whether the addition of HCN sig- 
nificantly improves the fit, we performed an F-test 
for an additional term (see lBevington fc Robinson 
l2003h . We find that for HCN, F^ = 0.993, and the 
probability of achieving this value for F^ with the 
addition of a random factor in our model is 32%. 
Thus, HCN does not significantly improve the fit 
at this wavelength range and the column density 
we find should therefore be considered an upper 
limit. Thus, the upper limit for HCN/HC3N w 1. 

We do not detect any absorption at the posi- 
tion of the bending mode of HC5N at 15.57 /im. 
Since the band strengths for the bendi ng modes 

Benilan et al. 2007 



of HC5N (268.2 cm-^atm-i 



and HC3N (245.1 cm^^atm'^ IJoUv et al.l 120071 ) 
are fairly similar and since the modes should have 
a similar profile, this suggests that HC5N is simply 
not present in SMP LMC 11. From the measured 
equivalent width of the band (0.007 /im) and the 
signal-to-noise ratio of ~85 in this part of the spec- 
trum, we estimate an upper limit of log N = 15.4 
to the column density for HC5N, yielding a lower 
limit to the HC3N/HC5N ratio of 10. 



5.4. Other Species 

Our best fit model shows a clear contribution 
from CH3C2H at a relatively high column density, 
blended with the C4H2 absorption band. Again, 
we performed an F-test and found that in this case, 
adding CH3C2H does indeed significantly improve 
the fit: We find F^ = 53 and the probability of 
observing this F-^ value with the addition of a ran- 
dom factor is ~10~^%. This absorption cannot be 
due to isotopologues of C4H2 either, as these ab- 
sorb at longer wavelengths than the main isotopo- 
logue peak (at 15.95 and 15.93 /tmforl£fCCCCH 



and HC^^cccjj respectively, I.Tollv et aLllioiol) 



We thus conclude that CH3C2H is indeed present 
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in the spectrum of SMP LMC 11. Note that this 
species was also observed in the pPN CRL 618 
at millimeter wavelengths, and that it was sug- 
gested to contrib ute to the absorption at infrared 
wavelengths too ( Cernicharo et al. 2001b[) . Here, 
however, we find a higher column density. 

Finally, we also find good fits and fairly high 
column densities for the CH4 absorption at 7.7 
lim and the C2H4 absorption at 10.53 /im. 

Judging from the residuals and low overall x^, 
we have accounted for most of the molecular ab- 
sorption. Thus, any additional molecular bands 
apart from those already noted must be either 
weak or perhaps shallow and broad. 

6. Discussion 

The rich molecular spectrum of SMP LMC 11 
seems to offer a unique astrophysical laboratory 
to study chemical pathways in carbon-rich envi- 
ronments, including the formation of benzene. In- 
deed, although the CSE of SMP LMC 11 shares 
some properties with CRL 618 — the only other 
pPN in which benzene is detected — there are some 
significant differences that offer clues to the con- 
ditions required for the efficient formation of ben- 
zene and other carbonaceous molecules. 

6.1. Circumstellar geometry 

The first aspect we consider is the geometry 
of the circumstellar environment and the physi- 
cal conditions of the material within. The con- 
tinuum emission in the Spitzer-IRS observations 
is due to dust and provides the background in- 
tensity against which the molecular gas absorbs. 
Thus, the molecular gas is either mixed in with 
the dust or located further from the star. Addi- 
tionally, the dust must be optically thick at in- 
frared wavelengths to explain the overall shape of 
the combined SL and LL spectrum. Moreover, op- 
tically thick dust is consistent with the featureless 
shape of the dust continuum which can be repre- 
sented reasonably well with blackbody curves (see 
also Section [3]). 

It is clear that the physical conditions in the 
CSE of SMP LMC 11 are somewhat different 
from those in CRL 618: the CSE of SMP LMC 
11 is denser and warmer than that of CRL 618. 
We find a typical dust temperature of 425 K, 
whereas the dust temperature of CRL 618 is 98- 



110 K (jPonfrfa et al.ll201ll) . Similarly, the molec- 
ular excitation temperatures are slightly higher — 
250-400 K versus 200-250 K— and we also find 
higher column den sities (e.g. N(C2H2) = 2 x 10^^ 



cm in CRL 618, ICernicharo et al.ll2001al versus 



> lOi* cm-2 in SMP LMC 11). 

If the CSEs of both CRL 618 and SMP LMC 
11 were simple outflows, the higher temperatures 
in SMP LMC 11 would indicate that the dust is 
located closer to the star and thus that SMP LMC 
11 turned off the AGB more recently. Indeed, 
the dust and molecular gas temperatures in SMP 
LMC 11 are more typical of a late AGB star than 
a pPN for which dust temperatures of 150-300 K 
are expected jKwok ,2000 ). However, as seen in 
Fig.H the Spitzer-IRS spectrum of SMP LMC 11 
also exhibits weak emission from [Ne II], which 
is typically one of the flrst excitation lines seen in 
young PNe. This clearly indicates that the central 
star is much hotter than a typical AGB star and 
is thus in the pPN stage. This apparent contra- 
diction is easily reconciled if we consider that the 
[Ne II] line does not originate from the same region 
as the molecular absorption and the dust. Indeed, 
since the dust is optically thick, the gas from which 
the [Ne II] originates cannot be located in a direct 
line of sight toward the central star: the dust hides 
the innermost region from view and no ionizing 
radiation can be expected to penetrate this thick 
dust layer. Thus, the [Ne II] line must originate 
from a region which is geometrically distinct from 
the location of the dust and the molecular gas and 
which is not obscured by dust. 

The current results suggest that the molecu- 
lar gas and dust surrounding SMP LMC 11 is lo- 
cated in a thick torus that we see fairly edge-on, as 
previously suggested in Paper I. In such a geome- 
try, the gas that causes the [Ne II] emission would 
then be in a bipolar outflow. The properties of 
the central object thus correspond to those of a 
pPN, but the circumstellar material does not cor- 
respond to a pPN outflow and evolves on its own 
timescale since it resides in a massive disk. Such 
a disk with optically thick dust would also shield 
the circumstellar material from the intense radi- 
ation of the central object. Note that a radially 
constrained dusty torus could also explain why a 
dust model with a constant-velocity outflow at a 
constant mass-loss rate overestimates the flux at 
the longer wavelengths. 
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Thus, from the geometrical point of view, SMP 
LMC 11 is very similar to CRL 618 for which a 
similar geometry involving a dense torus (in ad- 
dition to a bipola r outflow) is obseryation ally es- 
tablished (see e.g. Burton fc Geballe 19861 ). Such 
massive, vertically extended and long-lived disks 
are often seen around binary systems containing 
an evolved star such a s the Red Rectangle (see e.g. 
Jura et a"l] []J95: Waters et al. 111998), for example. 
However, it is not clear what is the primary cause 
for the higher temperatures in SMP LMC 11. This 
could simply indicate that the torus is closer to the 
central star than for CRL 618, either because the 
CSE is younger, or because it is expanding slower 
than that of CRL 618 (whi ch is expanding at a 
rate of 20 km s~^, see e.g. iHerpin fc Cernicharo 
20001) . However, one could also attribute the dif- 
ferent temperatures to differences in the optical 
thickness of the dust between the two objects. In 
either case, the molecular excitation seems con- 
strained to a fairly small range in temperatures. 

6.2. The chemistry in SMP LMC 11 

In terms of the molecular composition, the most 
striking difference between SMP LMC 11 and CRL 
618 is the much stronger absorption that we find 
for benzene. The center of the benzene band at 
14.85 fj,m is about three times deeper in SMP LMC 
11 than in CRL 618. It may be somewhat sur- 
prising then that we find a column density that 
is 13 times larger than wha t was found for CRL 
618 (jCernicharo et al.ll2001al ): however, we believe 
that this difference is due to the fact that we used 
more recent line lists. Additionally, not only is 
benzene itself much more abundant in SMP LMC 
11, it is also enhanced compared to most other 
species, especially compared to acetylene. If we 
take the lower limit for the column density of C2H2 
at face value, we find that there is about twice as 
much C2H2 as there is benzene; for CRL 618, the 
ratio found between these molecules was about 40 



(|Cernicharo et al.ll2001a^ . 

The relative abundances of the acetylene chains 
might offer some hints about the ch emistry lead- 
i ng to efficient benzene formation. IWoods et al 



(|2002l 12003^ fou nd that the interste llar benzene 
formation route ( McEwan et al. 19991) was not ef- 
ficient in environments such as CRL 618. Instead, 
they proposed a much more efficient route in such 
environments starting from acetylene and HCO^. 



The same model also predicts high column den- 
sities of C4H2 and C6H2 that are similar to one 
another. In fact, this matches the abundances in 
CRL 618 to within a factor 4-6. In the spectrum of 
SMP LMC 11 however, we cannot reliably estab- 
lish the presence of C6H2 in the spectrum. This 
suggests that there is yet another formation route 
to benzene in these warm and dense environments, 
or possibly some additional reactions that need 
to be considered. These reactions might involve 
fast and efficient ring-closing reactions that de- 
plete C6H2 or an additional route to benzene for- 
mation that starts from the abundant C2H2 and 
C4H2 which then decreases the efficiency of the 
formation of C6H2. 

It is equally interesting to consider the benzene 
loss reactions. Chemical models for environments 
comparable to those studied here show that the 
main destruction route for benzene is a further 
reaction with CN t o form bcnzonitrile (CgHsCN, 



Woods et al.l 120031 ). This molecular species has 



strong absorption bands at 13.2 and 14.4 /im, 
which we do not observe in the spectrum of SMP 
LMC 11. Although the CN in these reactions is 
expected to originate from the photodissociation 
of HCN which is again at best marginally present 
in the spectrum of SMP LMC 11. Benzene de- 
struction could thus be inhibited by the absence 
of the CN parent species. 

The cyanopolyynes al so present an interest- 
i ng ca se. As shown in ICherchneff fc Glassgold 
(|1993l ). for example, the primary routes to the for- 
mation of cyanopolyynes in pPNe environments 
are through reactions between members of the 
acetylene family and the CN radical which it- 
self is produced from HCN in photochemical re- 
actions. Reactions between CN and C2H2 then 
result in HC3N, while reactions with C4H2 yield 
HC5N. Since both C2H2 and C4H2 are abun- 
dant in the spectrum of SMP LMC 11, we would 
thus expect fairly large abundances of the longer 
cyanopolyynes as well. However, while we clearly 
detect a strong absorption band due to HC3N, we 
do not find much evidence for the parent molecule, 
HCN, nor for the longer cyanopolyynes. HC5N, 
for instance, should have a stron g bending mode 
at 15.57 /im (|Benilan et al.ll2007t ) which does not 
appear in the spectrum (see Fig. [5]). 

Such abundance patterns are very different 
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from those observed in CRL 618 and cannot be 
accommodated by the chemical models for these 
enviro nments (e.g. ICherchneff &: Glassgo ldl ll993l : 
Woods et alJl2003l) . Even if one would consider 
that HCN could be completely depleted by HC3N 
formation, one would expect to see efficient for- 
mation of HC5N as well since this involves the 
same mechanism and large amounts of C4H2 are 
available for this process. This observation then 
suggests that for the cyanopolyynes, some of the 
chemical pathways that are possible in pPNe en- 
vironments might be missing from the models. 

6.3. Pathways to PAHs? 

The formation of benzene is often studied in 
the context of PAHs, for which they are the ba- 
sic unit. In the formation pathways for PAHs, 
the formation of b enzene is considered to b e the 
bottleneck (see e.g. Allamandola et al. I ll989l) . We 
searched for any features due to neutral naphtha- 
lene and pyrene (two of the smallest PAHs consist- 
ing of two and three aromatic rings respectively) 
but did not find any evidence for these species. 
However, since PAH formation cannot occur at the 
low temperatures in SMP LMC 11 — typically 900- 
1100 K is required dFrenklach fc Feigel son 1989t 
Allamandola et al.l \l98^ ICherchneff et al.. ,19921) . 
this should not be surprising. 

It is interesting to consider what might hap- 
pen to the benzene as the central star and CSE 
evolve further. All other things being equal, the 
increasingly hotter central object would heat up 
the CSE. If the torus expands (as is the case for 
CRL 618), then the circumstellar material will be- 
come diluted and more transparent, which would 
clearly increase the importance of photochemistry. 
It is not immediately clear what the result of 
this increased photochemistry would be, but it 
is certainly possible that the conditions created 
are ideal for further processing of the circumstel- 
lar benzene into PAHs. If so, this might represent 
an important PAH formation pathway, provided 
that the torus contains enough mass to represent a 
significant carbon reservoir. Note that for objects 
like the Red Rectangle, the circumbinary disks are 
indeed found to be massive. However, in that par- 
ticular object , PAHs are clearly fo und in the bipo- 
lar outflows ( Waters et al. 1998[ ) and not in the 
torus. 

Although the molecular composition of the CSE 



of SMP LMC 11 is certainly unique thus far, we 
do not believe that the environment represented 
by this object is necessarily exceptional. CRL 
618 is the only other source where benzene is de- 
tected, and it also shows different abundances for 
the polyynes and the cyanopolyynes. However, 
the two environments are very similar in many 
of their physical properties. Given the short ex- 
pected timescales for the evolution from the tip of 
the AGB to the PN phase, it is reasonable to ex- 
pect the evolution of CSEs to be fast. The two 
objects could then represent cases where either 
the initial conditions (of the torus, for example) 
were slightly different or, alternatively, they could 
represent slightly different steps in the evolution 
of the torus. It would certainly be interesting to 
study how the physical conditions of a (slowly ex- 
panding) dense torus change in this short evolu- 
tionary phase and how the chemistry in this envi- 
ronment will evolve. 

7. Conclusions 

We presented an analysis of the rich molecular 
absorption in the Spitzer-IRS spectrum of SMP 
LMC 11. We have compared it to chemical mod- 
els for carbon-rich pPNe and to CRL 618, the only 
other pPN in which benzene is detected to date. 
The geometrical configuration in both objects is 
fairly similar and includes a dense, warm torus 
of material in which a rich molecular gas resides. 
However, the absolute and relative chemical abun- 
dances of the carbonaceous species in SMP LMC 
11 do not match models and are also quite differ- 
ent from those in CRL 618. In particular, ben- 
zene is very abundant, making SMP LMC 11 an 
important environment to consider in the study of 
benzene and possibly also in PAH formation. Cur- 
rent chemical models can reproduce some of the 
molecular absorption, but not all of it. In partic- 
ular, the absence of C6H2 and HC5N seems to re- 
quire additional ring-closing reactions that deplete 
those species. Alternatively, new chemical path- 
ways have to be found for the formation of C4H2 
and HC3N that do not result in abundant forma- 
tion of C6H2 and HC5N. We encourage a more de- 
tailed study of the chemical pathways which could 
explain these anomalies. 
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Benzene Bands at 6.8 and 14.8 fim with Model Fits 
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Fig. 3. — The benzene bands at 6.8 (top) and 14.8 
yCtm (bottom), the spectra are in black while the 
benzene models are in grey (red). The dotted- 
dashed vertical lines show the extent of the main 
absorption in the spectra, the dotted lines show 
the width of the model absorption for the entire 
band (at 6.8 /xm) and for the Q-branch (at 14.8 
lim). The dashed line in the top figure shows the 
center of the absorption, coincident with the emis- 
sion at the center of the 6.8 /zm band. See the 
electronic version of the Journal for a color figure. 



This 2-column preprint was prepared with the AAS lAT^TjX 
macros v5.2. 
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SMP LMC 1 1 and CRL 618 
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Fig. 4.— The normalized spectra of SMP LMC 11 
(black) and CRL 618 (dark grey or red in the elec- 
tronic edition) are shown with the C6H2 models. 
The model calculated using i nolecu lar constants 
from McNaughton fc Brugetl ( 1991 ) is shown in 
the solid light grey (blue in the electronic edi- 
tion) lines, while the modified model is in dashed 
Ught grey (blue) lines. The band center for the 
model calculated with molecular constants from 



McNaughton fc Brugetl (|l99ll ) is highlighted with 
a light grey (blue) dotted-dashed line. The dashed 
black line indicates an alternate fringing method 
for SMP LMC 11 wherein the appearance of the 
16.07 /xm feature is greatly diminished. See the 
electronic version of the Journal for a color figure. 
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